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Semiconductor quantum dots in silicon are promising qubits because of long spin coherence times 
and their potential for scalability. However, such qubits with complete electrical control and fidelities 
above the threshold for quantum error correction have not yet been achieved. We show theoretically 
that the threshold fidelity can be achieved with ac gate operation of the quantum dot hybrid qubit. 
Formed by three electrons in a double dot, this qubit is electrically controlled, does not require 
magnetic fields, and runs at GHz gate speeds. We analyze the decoherence caused by 1/f charge 
noise in this qubit, find parameters that minimize the charge noise dependence in the qubit frequency, 
and determine the optimal working points for ac gate operations that drive the detuning and tunnel 
coupling. 


I. INTRODUCTION 

Silicon semiconductor quantum dot qubits are promis¬ 
ing for quantum information processing because of their 
long spin coherence times and their potential for scala¬ 
bility and integration with classical electronics.® While 
high fidelity spin qubits in silicon have recently been 
achieved pEl for a practical quantum computer, it is also 
desirable to have a qubit that is purely electrically con¬ 
trolled, that does not require magnetic fields, and pos¬ 
sess fast gate speeds. The quantum dot hybrid qubitP® 
has the potential to meet these criteria. This qubit, 
which consists of three electrons in a double dot, has a 
qubit frequency of 11.5 GHz, set by the single dot singlet 
triplet splitting, and allows complete qubit control via 
the detuning-the voltage bias between the dots, and the 
tunnel coupling between the dots. Recent experiments 
have demonstrated ^ 85% gate fidelities for dc (direct 
current) operation and ^ 93% for resonant ac (alternat¬ 
ing current) operation that drives the detuning. The gate 
operations that limit the qubit fidelity in the these exper¬ 
iments involve transitions between qubit states, which we 
call X rotations. While ac gate operations, which drive 
these transition in a manner similar to electron spin reso¬ 
nance, have resulted in a significant improvement in qubit 
fidelity, fidelities above the threshold for quantum error 
correction have not yet been achieved. 

To understand the limiting factors in ac gate oper¬ 
ations, one needs to examine decoherence for driven 
qubits, which is substantially modified from the case of 
free qubit evolution. This subject has been studied ex¬ 
tensively both theoretically and experimentally for su¬ 
perconducting qubits P® l n contrast to free evolution, 
where typically all noise power below the precession fre¬ 
quency contribute to dephasing, during driven evolution, 
the effect of low frequency noise on the qubit dynamics is 
mitigated, resulting in significantly improved coherence 
times. On the other hand, the driven qubit is exposed to 
noise at specific high frequency components in the noise 
spectrum. In particular, dephasing of ac resonant gates 
for the hybrid quantum dot qubit is sensitive to the noise 
power at the qubit frequency. 

In this paper, we perform a systematic analysis and 


numerical optimization of ac resonant gate fidelities for 
the quantum dot hybrid qubit. We develop a decoher¬ 
ence model for the hybrid qubit that fully takes into ac¬ 
count the 1/f charge noise spectrum, which is the domi¬ 
nant source of decoherence in double quantum dots. Our 
model also takes into account nonlinear qubit dynamics 
that occur under strong driving conditions. While we 
consider only the hybrid qubit in this paper, our model 
can be readily applied to other semiconductor qubits. 

Our approach to optimization has two parts. First, 
we consider dephasing effects due to the detuning de¬ 
pendence of the qubit excitation frequency, which we 
call charge dispersion , borrowing terminology from the 
superconducting qubit literature. Charge dispersion de¬ 
scribes the sensitivity of the qubit frequency to charge 
noise fluctuations in the detuning and associated dephas¬ 
ing. I n a sp irit similar to the transmon superconducting 
qubit J® 12 we reduce the charge dispersion of the quan¬ 
tum dot hybrid qubit by tuning the static tunnel cou¬ 
plings. Recent experiments with parameters approaching 
this optimal parameter regime has improved X-rotation 
coherence times from 33 ns to 177 ns. 13 

However, because detuning is both the dominant noise 
source®! 17 and a drive parameter in the hybrid qubit, 
minimizing charge dispersion while maintaining fast gate 
speeds, necessary for high fidelity operation, is problem¬ 
atic. We show that this problem can be circumvented by 
driving instead the tunnel coupling, which is more effi¬ 
cient, and results in higher gate speeds and fidelities than 
driving the detuning. 

Second, we numerically optimize the X-rotation fi¬ 
delity as a function of the ac drive amplitude and de¬ 
tuning. In this optimization, we simulate qubit dynam¬ 
ics using an effective two-dimensional Hamiltonian, to 
which we apply the Bloch-Redfleld master equation for 
driven qubits and analytic formulas for dephasing rates 
due to 1/f noiseP® 18 19 We consider both the case of 
detuning and coupling drive, and find that X-gate fideli¬ 
ties exceeding 99% can be achieved in both cases, and 
fidelities of 99.8% can be achieved by driving the tunnel 
couplings. This optimal fidelity agrees with the result of 
simulations using a three state Hamiltonian that includes 
the nearest leakage state, averaged over numerically gen- 
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FIG. 1. (Color online) Illustration of the ac driven hybrid 
double quantum dot qubit. The spin configuration of the 
state | -S) and the single dot energy levels corresponding to 
singlet and triplet state occupation on the right dot are shown. 
The electrostatic confinement potential is characterized by the 
detuning e 0 , the potential energy difference between the left 
and right wells, and the barrier height, which sets the tunnel 
coupling (Ai,A 2 ). Detuning and tunnel coupling drive are 
implemented by modulation of the detuning (e ac ) and the 
barrier height, respectively, shown in dashed lines. 



FIG. 2. (Color online) Qubit energy levels Eq and E\ and 
the leakage energy level E 2 as a function of detuning e com¬ 
puted with the three state Hamiltonian in Eq. (TlJ) for param¬ 
eters (Est, Ai, A 2 ) = (50,10.8,30.6) fieV taken from Ref. [21] 
Black dashed lines: qubit energ y lev els computed from the 
static effective Hamiltonian Eq. (A10) 


erated 1// detuning noise. 

This paper is organized as follows. In section [TlJ we 
summarize the relevant features of the quantum dot hy¬ 
brid qubit. In section m we explain our approach to 
finding the optimal static tunnel couplings that minimize 
charge dispersion. In IV and [V| we describe our model¬ 
ing of ac resonant gates and decoherence, respectively. In 
section |vg we present the results of our numerical sim¬ 
ulations on ac gate fidelities as a function of detuning 
and drive amplitudes. A general procedure for deriving 
the effective Hamiltonian for a driven qubit is given in 
appendix [Aj and a summary of pure dephasing and re¬ 
laxation rates due to 1// noise, as well as a numerical 
procedure for numerically generating such noise, is given 
appendix [C] 


II. QUANTUM DOT HYBRID QUBIT 

The “quantum dot hybrid qubit”®^ is formed by 
the manifold of three-electron states in a double quan¬ 
tum dot with total spin quantum number 5 = 1/2 
and ^-projection S z = —1/2. The three-electron dou¬ 
ble dot Hubbard Hamiltonian is given in Ref. 0J 
A basis for the relevant states in the hybrid qubit 
regime can be chosen as |-5) = |^5), |5-) = |5^), and 
|-T) = y/l/3\iT 0 ) + v/273|tT_), where S(T) refers to 
singlet (triplet) states on one dot. In our nota- 
tion, | IS) s= denotes spin down electron on 

the left dot and a singlet on the right dot, and 
|t T~) = Jt) L | T-) R , etc. The spin states of |-5) and |-T) 
are identical to those of the exchange-only logical qubit 
which consists of three electrons in a triple dot, instead 
of a double dot. This qubit is operated typically at an 
electron temperature of T — 140 mK. 21 

Fig. 0 illustrates the basis states {|-5), |-T), |5-)} in a 
gate-defined electrostatic potential. The control parame¬ 


ters of the quantum dot hybrid qubit are the double dot 
detuning e, defined as the energy difference between (2,1) 
and (1, 2) charge states, and the tunnel coupling Ai (A 2 ) 
which cause charge transitions between |5-) and |-5) (|5-) 
and |-T)). The Hamiltonian in the basis {|-5), |-T), |5*)} 
is given by 


e (° 0 M 

EL(e, Ai, A 2 ) = —- + 0 Est —A 2 1 , (1) 

2 \A X -A 2 6 J 

where Est is the energy splitting between the lowest ly¬ 
ing singlet and triplet states on the right dot. The en¬ 
ergy levels of the Hamiltonian Eq. 0 , denoted by | 0 ), 
| 1 ) and 12 ), ordered from low to high energy, are plot¬ 
ted in Fig. [2] as function of detuning, for the parameters 
(Est, Ai, A 2 ) = (50,10.8,30.6) /ieV taken from Ref. [ 21 ] 
The qubit logical states are |0) and 11), while |2) is a 
leakage state. We will represent the qubit on a Bloch 
sphere where | Z) = |1) and | —Z) = |0) are located on the 
north and south pole, respectively. In this work, we con¬ 
sider operations of this qubit in the ( 1 , 2 ) charge regime, 
e > Ai and e — Est > A 2 , where it is useful to think of 
the qubits states as mainly comprising of singlet-triplet 
states on the right dot, |0) ~ |*5) and |1) |-T), with a 

small hybridization of the |5-) state. 

The tunneling couplings Ai, A 2 g enerally have expo¬ 
nential dependence on detuning 4 22 4 In order to fit the 
experimentally observed resonant frequencies in Ref. [ 2 j.j 
we find it necessary to introduce such a dependence into 
the second tunnel coupling as 

A 2 -» A^(e) = A°e“ |e “ BsT|/e2 , 

with 62 = 400 peV, while A° has a sufficiently long expo- 
nential decay length that it can be regarded a constant 
independent of e. This detuning dependence is included 
in Fig. [2] and in all results reported in this paper. 
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FIG. 3. (Color online) As function of detuning e: (a) qubit 
energy levels Eq and E±, (b) qubit frequency cjz/27t = 
Eio/h, and (c) the quasistatic dephasing time scale T 2 *, for 
(A?,A°) = (10,30) fieV (Dashed, blue) and (A?,A°) = 
(30,30) peV (Solid, Maroon). 


III. OPTIMAL STATIC TUNNEL COUPLINGS 

In this section, we tune the static tunnel coupling to 
minimize charge dispersion, defined as the detuning de¬ 
pendence of the qubit frequency cjz/27t = Em/h , where 
Eio = Ei — Eq. In the first instance, this will suppress 
the pure dephasing rate of free qubit precession (Z rota¬ 
tions) due to quasistatic detuning noise, given to leading 
order by y/2/T^Z) = cr* | duJz/de |, whertP^ 


a * = c ^l ln & =5 - 7fieY (2) 


is an effective quasistatic noise variance, c e = 2.38 /ieV is 
the parameter in the 1// detuning noise spectral density 
S e (uj) = c 2 /uj [c.f. Eq. (19)], determined from experi¬ 
mental data on T\ relaxation, as described in appendix 
|C 2[ and uji « 1 Hz is the low frequency cutoff for the 
1// noise. This cr* value is consistent with the static 
noise variance previously used to model the experiment of 
Ref. |7j A detail discussion of pure dephasing rates due to 
quasistatic noise, including logarithmic correctio ns c har- 
acteristic of the 1// noise, is given in appendix C3 and 
Eq. § is derived in appendix |C 3 b| 


At the same time, minimizing charge dispersion im¬ 
proves the coherence time of ac-driven X rotation, which 
also suffers quasistatic dephasing from low frequency fluc¬ 
tuations of the Rabi frequency due to charge dispersion at 


quadratic order where Sujz = SeduJz/de , Se denotes 
detuning noise, as discussed in section [YD) Z rotations 
are also important because they can be combined with 
ac driven X rotations to achieve universal single qubit 
control. 

As evident from the energy level diagram, deep in the 
(1,2) charge regime (e Ai, A 2 ), the qubit frequency is 
insensitive to the detuning. Charge dispersion is minimal 
in this regime because here logical states differs mainly 
in their spin instead of charge character. As a result, 
free precession is protected from dephasing due to charge 
noise fluctuations in the detuning. A recent experiment 
in this regime demonstrated free induction decay at a 
frequency of 11.5 GHz with a dephasing time of 10 ns, 
resulting in a Z-gate fidelity of 96%We next show that 
this fidelity can be significantly improved by tuning the 
static tunnel couplings. 

To find the optimal static tunnel couplings parame¬ 
ter regime, consider the following qualitative argument. 
Charge dispersion comes from level repulsion at the an- 
ticrossings due to the tunnel couplings. Specifically, as 
one moves towards e = 0 from large detunings e Ai, 
the energy level E\ is repelled downwards at e = E$t by 
A 2 and then repelled upwards at e = 0 by Ai, see Fig.[2j 
We thus expect the net effect of level repulsion to be 
minimized when Aj ~ A®. By performing an empirical 
search, we find minimal charge dispersion at the tunnel 
couplings A? = A° = 30 /ieV. Fig. [ 3 ] a-c compares the 
qubit energy levels, excitation frequency Co>z/27r, and the 
quasistatic dephasing time scale T| (Z) for the optimal 
tunnel couplings with the ones reported in Ref. [24j The 
curvature in the energy levels near the anticrossing and 
associated charge dispersion is greatly reduced, and, as a 
result, T 9 * (Z) is increased by an order of magnitude from 
10 ns to 180 ns. 

In principle, the high frequency components of the 
1/f noise spectrum can also cause dephasing, due to 
quadratic charge dispersion 5e 2 d 2 uJz /<9e 2 P 9 The result¬ 
ing decay envelop is exponential with a decay time r 2 
inversely proportional to the curvature of the charge dis¬ 
persion t 2 (Z) oc (d 2 ioz/de 2 )- 1 . However, this decay is 
strongly suppressed for the optimal tunnel couplings, as 
t ~2 > 1 /is in the far detuned regime, as shown in appendix 
|C3c| Since the the Z-gate times are about t g (Z) ~ 0.1 
ns, this decay should lead to a Z-gate infidelities less than 
t g (Z) /t 2 (Z) ~ 0.01%, which is very small. 

Henceforth, we set static tunnel couplings at their op¬ 
timal values, and optimize ac X gates as a function of 
remaining parameters. 


IV. AC RESONANT GATES 

Ac resonant techniques for producing transitions be¬ 
tween the logical qubit states (X rotations) are similar 
to electron spin resonance (ESR), where the logical qubit 
plays the role of the electron spin and oscillations at the 
qubit frequency ujz / 27r in the electrostatic control param- 
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eters plays the role of the resonant driving field. Even in 
the absence of decoherence, several complications arise 
in “logical qubit resonance” (LQR)^ which do not occur 
in ESR. In implementing LQR, one has indirect control 
over the direction and magnitude of the driving field, 
which is determined by the qubit response to perturba¬ 
tions in the control parameters. Furthermore, to achieve 
fast gate speeds, one may have to enter the strong driv¬ 
ing regime where the qubit’s nonlinear response can spoil 
gate fidelity. While this problem could be addressed by 
properly shaping the driving pulse, as demonstrated with 
diamond NV centers pMSH we no t consider these tec- 
niques in this paper. 

LQR for the hybrid qubit is performed by adiabatically 
initializing to a detuning eo, and then applying ac oscilla¬ 
tions of the detuning (voltage bias between the quantum 
dots) 


e(i) = e 0 + tacit ), 

Ai(i) = A? 

A 2 (t) = A° 2 e~^ t) ~ EsT)/e2 , (3) 

or tunnel couplings (DQD barrier height), 

Ai(t) = # A ac (t) 

A 2 (£) = A 2 (eo) + A ac (t) , (4) 

where we will consider the specific drive signals, 

(aTc ( S ) =cosuzt (£) ’ =E ° i/h • 

where A e and Aa are the drive amplitudes of detuning 
and coupling drive, respectively. These two types of driv¬ 
ing are illustrated in Fig. [I] Note that, due to the energy 
dependent tunnelling, detuning drive will also result in 
an ac signal in the second tunnel coupling, as well. The 
total Hamiltonian now comprises of a dc and ac compo- 
nent, Hit) = Hdc + W-adt), where U dc = %(e 0 , A^A-j) 
is the dc component, and 

n ac (t)=H[e(t),A 1 : 2 (t)]-n dc 

is the ac component. 

Tunnel coupling driving should enable faster, more effi¬ 
cient (fast speed per drive amplitude), and higher fidelity 
gates than detuning drive, because of stronger transition 
matrix elements between qubit logical states at the large 
detunings where the qubit is protected from dephasing. 
There, charge hybridization between (1,2) and (2,1) are 
minimal, so detuning drive, which do not couple definite 
charge states, cannot be effective, while transitions be¬ 
tween |0) — |-5) and |1) ~ |*T), mediated by tunneling 
(see Fig. [l]), have matrix elements of the order AiA 2 /e. 
This drawback of detuning drive is quite generic: since 
the gate speed scales as the double dot susceptibility to 
detuning fluctuations, increasing drive efficiency also in¬ 
creases sensitivity to detuning noise. This is particularly 


important for finite frequency noise, because the rotating 
frame pure dephasing rate 1/T£ for ac X rotation scales 
quadratically with the detuning drive Rabi frequency, as 
shown in section V C, Eq. ( [23]) . This problem is evi¬ 
dent in the experiment of Ref. 121 [ where high speed (^ 1 
GHz) ac gates were achieved at detunings near the (1, 2) 
to (2,1) charge transition, but the decay times were only 
1-2 ns, even at the charge noise sweet spot, limiting gate 
fidelities to ~ 86%. 

To model ac gates, it is more convenient to ex¬ 
press the Hamiltonian in the basis of energy eigen¬ 
states, ii(eo) = (eo)l-LU(eo), where U. is the unitary 
transformation between eigenstates {| n)} of the static 
Hamiltonian at eo and {\-S), |-T), \S-)}. Using this 
Hamiltonian, we can estimate the leakage probabilities 
for short times from qubit states | m) to the leakage 
state |2) by p rn ^2 — (; Hm 2 /E 2 m ) 2 , where m = 0,1 and 
£’ 2m = F 2 — E m are the energy gaps to the leakage state. 
The leakage is very small at large detunings, < 0.1% for 
e > 200 /ieV, due to the large energy gap to the leakage 
state, £’ 2m > 150 peV. We have also verified this leak¬ 
age estimate by numerical solution of the density matrix 
equation using the three-state Hamiltonian Eq. 0- Since 
the leakage is negligible, the qubit dynamics is governed 
by an effective two dimensional Hamiltonian, presented 
in section HV Ai 


A. Effective Hamiltonian for the driven qubit 


In this section, we analyze the qubit dynamics and 
noise using the effective Hamiltonian for the driven qubit, 
which is useful for both numerical and analytical calcula¬ 
tions, for providing intuition, and for applying techniques 
of electron spin resonance to the driven qubit. Applying 
standard techniques 29 30 described in appendix[Aj we find 
an effective Hamiltonian given by a perturbative expan¬ 
sion in the inverse energy gaps between the qubit and 
the leakage state, E^i and E% q 1 , and can generally be 
written as 


h = b(e 0 , e ac (t) + 6e(t), Ai, 2 (t)) • a/2 . 


( 5 ) 


where Se(t) is the detuning noise, a are the Pauli matri¬ 
ces, b is the effective magnetic field that is a nonlinear 
function of the qubit control parameters (e, Ai, A 2 ). We 
retain the leading ord er term in the effective Hamilto¬ 
nian, given in Eq. (A8). The static effective Hamiltonian 
given in Eq. (A10), was derived in Ref. 131, and the 
qubit energy levels computed with it, shown in Fig. [2j 
agrees very well with the ones calculated with the static 
three-state Hamiltonian Hdc [Eq. ([!])]. However, h gener¬ 
ally has a different functional dependence on eo and e ac , 
as indicated in Eq. 0, due to eo-dependent transforma¬ 
tions used to approximately block diagonalizes the static 
Hamiltonian ( Hdc ) between qubit and leakage states. 
These transformations define the basis in which hd c is 
derived, which i s rela ted to the basis {|-5), |-T), | S-)} of 
Eq. 0 by Eq. ( |A11| ). In particular, we emphasize that 
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one cannot derive Eq. © by simply taking eo —► eo + e ac 
in h dc . 

In the far detuned limit (e —> oo), the matrix elements 
of Eq. © are given by 


h _ e A i 

h°°, 


hoi — 


AiA 2 


hn = Est ~ - 


1 


Est — e 


Aj 

Est — £ 


( 6 ) 


The quadratic terms describe virtual transitions to the 
leakage state, which would otherwise appear as second 
order terms in the time evolution operator for the 3D 
Hamiltonian Eq. Q i l32 l 33 l 

The effective magnetic field Eq. © has three compo¬ 
nents 


b (t) = b 0 + h ac (t) + db(f), (7) 

as illustrated in Fig. [4^,, where b 0 = b(eo) is the static, 
dc field corresponding to Hdc , which will define the “lon¬ 
gitudinal” direction in the lab frame, b ac is the ac field 
corresponding to H ac , and 5b is the effective field due to 
noise, defined by b(eo + 5e) — b(eo). We can generally 
parametrize the static field as 


bo = sin + cos #oz , 


where the detuning dependence of the angle Oo = 0(e o), 
along with the relative angles between the drive field B ac 
and bo, for both types of driving, denoted by 0^} and 
0a^\ are plotted in Fig.^. Unlike conventional spin res¬ 
onance, b ac has both longitudinal and transverse compo¬ 
nent, but only the transverse field b±(t) = bo x b(f) x bo 
drives qubit transitions. 

We next present a general analysis based an expansion 
of the effective Hamiltonian up to second order in drive 
and noise fields that captures all the relevant physics in 
the parameter regime of interest. The results reported 


below in section VI, however, follow from numerical sim¬ 
ulations with the effective Hamiltonian Eq. (A8) that in¬ 
clude higher order nonlinear effects. 

The noise field due to fluctuations in the detuning e 
and the singlet-triplet energy splitting Est is given by, 


5b~5e 


db 

d(Se) 


5 Est 


db 

DEst 


(8) 


Note that this linear expansion will lead to second order 
noise correlations and associated relaxation terms in the 
density matrix master equation, see section [Vj which is 
derived by solving for the system-bath interaction to sec¬ 
ond order in perturbation theory? 34 The driving field up 
to quadratic order is given by 


bac = b^ + b£) 

biV = + ^b ac , (5b, 

h (2) = ^ 2 0) d 2 b 
2 du 2 


u(t)5e 


d 2 b 

dud(5e) 


(9) 


( 10 ) 



e (peV) 


C Lab Frame (LF), RWA d Rotating frame (RF), RWA 



FIG. 4. (Color online) (a) At a fixed detuning eo, X — Z axes 
of the Bloch sphere, representing the qubit logical basis, and 
the x — z axes, representing the basis in which the effective 
Hamiltonian Eq. © is derived. The Z axes lies along the 
static field bo, at an angle 0q rotated from the z axis. The 
ac drive field b ac lies at an angle 0 ac from the Z axis, and 
has the transverse component B ac x (blue arrow). 5b is the 
noise field with transverse component 5Bx, shown in green, 
(b) Detuning dependence of the static field angle Oo and the 


ac drive field angles for detuning (0^1) and tunnel coupling 
(eiS) drive, (c and d) Effective fields in the rotating wave 
approximation when a phase 0 is included in the drive ampli¬ 
tude, which goes as cos(c ozt + 0). (c) Lab frame: the trans¬ 
verse resonant component of the ac drive field B^ c (blue ar¬ 
row) and the noise in its drive amplitude 5~Bac (green arrow) 
[Eq. ©], is rotated at an angle uozt + (j) from the X axis at 
time t. Th e X and Z components of the noise field, 5Bx and 
H5ujz [Eq. (12)], respectively, are shown in green, (d) Rotat¬ 
ing frame: The transverse ac drive field is fixed at an angle 0 
and has magnitude Hujr (blue arrow) with a noise component 
H5uj r (green arrow). The tranverse noise field ^B ; ± (green 
arrow) [Eq. ©] counter-rotates at the qubit frequency ljz- 


where u = e ac or A ac , and biV and bi? are the term 
linear and quadratic in u. In the second equation, the 
term proportional to 5e represent fluctuations in the drive 
amplitudes due to detuning noise. 

In the Bloch sphere representation of Eq. ©, the qubit 
logical states, defined by the eigenstates of bo • cr, lies 
along the axis defined by bo. It will be convenient to use 
instead a basis in which the logical states always lie along 
the Z-axis. To this end, we apply the unitary transfor¬ 
mation 


_/cos(6» 0 /2) -sin(6» 0 /2)\ 

0 L in (V2) cos(# 0 /2) ) ’ 

which diagonalizes the static Hamiltonian, 
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FIG. 5. (Color online) Rabi frequency as a function of detun¬ 
ing e, computed from Eq. (13), for detuning drive amplitude 
A e = 30 jueV and coupling drive amplitude A a = 15 jueV. 


Uq (bo • cr)Uo = az- We write the effective Hamil¬ 
tonian in the local basis as 

H = U\e)hU{9) = \{B z (t)az + B x (t)a x ) , (11) 


where Bi(t) = b (t) • e$, {e^} are detuning-dependent 
qubit pseudospin basis vectors given by ez(co) = b 0 , 
e x(^o) = cos <9 0 x — sin#oz, and ey = y, and indices in 
capital letters denote local frame axes. These local axes 
are illustrated in Fig. [4^i. 

Up to th e exp ansions Eq. (10), the effective field com¬ 
ponents in (|A8|) are given by 


B z (t) = -h(u z + Su> z ) + B^ z + B™ z 

B x (t) = 2h{ui R + 5u R ) cos u> z t + SB X + B^J X . (12) 


where the Rabi angular frequency ujr and its fluctuation 
5ujr are given by 


Ur 


A u dBx 

2 h du 


5ujr = 5e 


du R 
dSe ' 


(13) 


The noise term Sujr is a nonlinear effect that comes from 
second order processes in which e ac drives a transition 
from the qubit subspace into the leakage state, and then 
noise Se drives transition back into the qubit subspace, 
or vice versa, as illustrated in Fig. Oh 

The Rabi frequency \ujr\/2tt is plotted for detun¬ 
ing drive with A e = 30 /ieV and coupling drive with 
A a = 15 /ieV in Fig. [5] As expected from our argument 
in section |IV[ tunnel coupling driving is more efficient 
and stronger at large detunings. On the other hand, fast 
detuning drive gate speeds can be attained at lower de¬ 
tunings where charge dispersion is larger. We can quan¬ 
tify this relationship by considering the scaling of the 
charge dispersion curvature with the detuning drive effi¬ 
ciency d Bx / d e ac = 2hc j 0 r/A 6 . From Eq. (13), Eq. (C5), 
and Eq. (Cl 1), we find 


d 2 E w 2H(co R /A e ) 2 

8e 2 K u z 


(14) 


where we used db/d(5e) = 3b/<9e ac P^ see Eq. ([5]), and 
neglected terms that are independent of ujr which comes 
from d 2 Bz/deQ and the difference between db/deo and 
db/de ac , see Eq. (A9) in appendix |A| The charge disper¬ 
sion thus scales quadratically with ujr. 


B. Rotating wave approximation 


While we will simulate the unitary qubit dynamics gov¬ 
erned by the full effective Hamiltonian in section |VI[ it 
will be instructive to analyze the qubit dynamics in the 
rotating wave approximation (RWA). The leading correc¬ 
tions to the RWA are of order ur/ujz, see, e.g. , Ref. HU 
The rotating frame (RF) Hamiltonian, defined by 

h rf = -ihuld t u z + U Z HU Z (15) 

where Uz=e ' luJztaz ^ 2 , is given in the RWA by 

HT A) = \ + Su r )<j x + 5u z a z \ + <5B(,_ • ^ 

+ (16) 

where the transverse noise field has additional time de¬ 
pendence (indicated by ') corresponding to Z rotations 
at the driving frequency uoz 

5B(j_(t) = SBx(t )(cos ujzt^x ~~ sincj^tey), (17) 


while 5(Jz and Sujr are unaffected by the transformation. 
Note that the addition of a phase p in the ac drive signal 
(coscj^t —»> cos (uzt+ (/))) enables rotations in the RF 
about an axis in the XY plane at an angle p rotated 
from the X axis, as illustrated in Fig. [4]i. Since only two 
axis control of the qubit is required for universal gate 
operations, this qubit could be operated entirely with ac 
gates. Fig. [4ji illustrates the relevant effective fields in 
the rotating frame. 

In addition to the standard spin resonance Hamilto¬ 
nian in the rotating frame, Eq. (16) has a longitudinal 

drive B ac z and a nonlinear drive term B K ac J . The longi¬ 
tudinal driving term does not affect our results because 
we operate the qubit in the regime ujr ~ uoz^ Details 
about the magnitude and effects of this term are given in 
appendix [B] The second harmonic of the driving field is 
given by 


r ( 2 ) = .2 1 + cos2w z £ <9 2 B 

ac u 4 q u 2 • 

The ac term with frequency 2 ujz is off-resonant, so it 
can be neglected consistently in the RWA. On the other 
hand, the dc term drives rotations at a frequency set 
by A^(d 2 b/du 2 )/4, which can cause beating with oscil¬ 
lations driven by the linear transverse ac driving field 
B ac x • This beating pattern can be seen in oscillations 
of the infidelities as a function of e for large drive ampli¬ 
tudes, as seen in Fig. of section [V C| 
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V. DECOHERENCE 


In the following, we apply the decoherence model for 
the driven qubit in the presence of If f nois e, previously 
developed for superconducting qubitiP"^^, to the quan¬ 
tum dot hybrid qubit. In section [V A[ we discuss the rel¬ 
evant noise sources and their power spectrum. In section 
VB[ we show how ac driving avoids the low frequency 
noise in SBx • In V C[ we apply the Bloch-Redfield equa¬ 
tions in the rotating frame to the driven qubit, and com¬ 
pute the associated relaxation rates. In|VDl we estimate 
the effect of the low frequency noise in the Rabi frequency 
Sujr. The results derived from our decoherence model are 
checked with simulations using the three state Hamilto¬ 
nian with numerically generated Iff detuning noise in 
section ED 


A. Noise sources, power spectrum, and relaxation 
tensors 

As we mentioned in the introduction, charge noise with 
Iff type noise spectrum is the dominant cause of de¬ 
coherence for double quantum dots in Si. Electrostatic 
coupling to this charge noise causes fluctuations in the 
detuning Se and the singlet triplet splitting SEst , which 
are completely characterized by the classical noise auto¬ 
correlation 


S e (t — t') = (Se(t)Se(t')) 

Ssr{t — t') = ( SEsT{t)SEsT{t‘')) • (18) 


in e and Est- They are plotted as a function of detuning 
in Fig. [6] a and b. The decrease in Gij(e ) as a function of 
detuning is consistent with our discussion in secion II and 
III: when the logical states have the same charge char¬ 
acter, Se cannot cause changes in the qubit frequency or 
transitions between logical states. On the other hand, as 
the qubit logical states acquire different spin character as 
e increases, they become more sensitive to SEst noise, 
so that Fij increases. 


We will next make some approximations to s lj ap¬ 
propriate for the optimal parameters considered in this 
work. First, detuning noise is by far the dominant 
noise parameter, with c e = 2.38 //eV, while an estimate 
in Ref. 14 T1 gives cst ~ 10 -3 /xeV. The large difference in 
noise strength occurs because the detuning couples to 
charge noise via the double dot electric dipole moment, 
but while the single-triplet splitting Est couples via a 
single dot quadrapole moment. 41 ! Since even the largest 
relaxation rate associated SEst , set by the noise power 
at the ujz or ujr as discussed in section |V C[ is given by 
FzzCst/^ur ^ 100 Hz, negligibly small compared to 
other relaxation rates considered in this work, we can 
neglect SEst altogether in the decoherence of ac gates 
Second, we can neglect the anisotropic relaxation term 
Gxz which represent noise correlations in different direc¬ 
tions on the Bloch sphere, as it is an order of magnitude 
smaller than the diagonal term Gxx as shown in Fig.[6^. 
Note that Gzz and Gxz are strongly suppressed because 
we already tune the parameters to minimize charge dis¬ 
persion ( dBz/de ). 


Their noise power spectrum S(u) = f dte lujt S(t) are 
given bjEl 


S e (w) = ^r, Sst(u) = ui<\ui\<ui h , (19) 
M M 


where we impose a sharp high and low frequency cutoff, 
u)h and ujr respectively. The low frequency cutoff uji is 
set by the total measurement time uji = 2tt /r m ^^, which 
is r m = 200 ms in the present experiment, corresponding 
to a low frequency cutoff of r” 1 =5 HzP^ 

Dephasing and relaxation of qubit dynamics, described 
by the master equation in section |V C[ is not determined 
directly by the noise power in Eq. (19), but the power 


spectrum of the noise field that appears in the effective 
Hamiltonian Eq. 0, given bjpsi 

Sij(u) = (5Bi(u})5Bj(—u})) = GijS e (io) + Fj^SsTiT) 

( 20 ) 


where we define the relaxation tensors 


_ dB t d_Bj OB, dBj 

lj[ ’ d(6e)d(Se)’ vU dE ST dE ST ' 


( 21 ) 


These tensors, which are anisotropic and detuning de¬ 
pendent, describe the qubit susceptibility to fluctuations 


While phonons cannot cause transitions between dif¬ 
ferent spins, they could potentially cause decoherence 
as they couple to the difference in the charge distri¬ 
butions between singlet and triplet states in the dou¬ 
bly occupied dot with an interaction term that goes as 
H ep ~ (a ph + «ph)(l' 5 ') (‘ s \ ~ l T ) (‘ T l)> where a ph is the 
phonon destruction operator.^ However, averaging the 
fluctuations of the relative phase factor between | -S) and 
| -T) [see Eq. ( |C14| )] over an incoherent thermal bath does 
not cause decay because the phonon density of states van¬ 
ishes at low frequencies,^ in contrast to the pure dephas¬ 
ing by 1// detuning noise. On the other hand, according 
to Ref. [43j when the phonon dissipative dynamics and 
coherent interaction with the qubit is fully taken into ac¬ 
count, phonon relaxation can cause exponential decay of 
the singlet triplet coherence. To model this effect, we 
will take the decay time T p h — 1 /is calculated in Ref. |4T| 
for silicon with a lateral electron confinement length of 
40 nm. Since the phonon coupling has the same form as 
the charge noise coupling to Est , the phonon relaxation 
tensor is Fij. In particular, for the detunings of interest, 
only Fzz — 1 is relevant. Therefore, in the following, we 
will directly incorporate the decay due to the phonon de¬ 
phasing in the off-diagonal density matrix elements by 
taking p 0 i(t) -> / o 0 i(i)e _t/Tph . 

























verse field 



e(/veV) 


FIG. 6. (Color online) For (A?, A°) = (30,30) /zeV and as 
a function of detuning e: (a and b) The relaxation tensors 
Eq. ( |2T| (a) Gij for detuning noise and (b) Fij for noise 
in the singlet-triplet splitting Est■ (c) For the drive ampli¬ 
tudes ( A e ,A g ) = (30,10) /xeV, relaxation times (T£, Tip,Tip) 
Eq. (26) in the Bloch equations for both types of ac drive, 
and the phonon relaxation time T p h- (d) Illustration of the 
longitudinal and transverse relaxation process in the rotating 
frame, relative to the drive field (ur), governed by Tip, and 
Ti P , respectively, and the phonon induced relaxation (T p h). 


B. Mitigating dephasing from low frequency noise 


A key advantage of ac gate operations is that it is in¬ 
sensitive to low frequency transverse noise in SBx , re¬ 
sulting in improved coherence during driven evolution 
compared to free induction decay. Quantitative analysis 
of this effect is more convenient in the rotating frame, 
where the qubit simply undergoes free rotation about 
the X axis, see Fig. App lying the pure dephasing 
rate formula given in Eq. (C8), and noting that due to 
Eq. ( fl7l ), the transverse noise ' ± {t) has a shifted noise 
spectrum given by — 


of r \ of ( \ Sxx(u + luz) + Sxx(u 
bxxyu) — &yy\ u) — ^ 

one finds the pure dephasing rate, 

1 _ _ Sxx{wz) 

T l ~ 2 h 2 ~ ' 


wz) 

1 

( 22 ) 

(23) 


Thus, dephasing due to the low frequency transverse 
noise in SBx is avoided. The physical origin for this 
effect can be understood by considering the RWA in the 
lab frame as illustrated in Fig. if where the qubit pseu¬ 
dospin precesses about the instantaneous, rotating trans¬ 


B_l (t) = (h(jjR cos oozt + SBx (t))ey + TujJr sin (jJzte Y • 

The relative phase S(j) accumulated due to fluctuations of 
the instantaneous rotation frequency is given by 

S(/)= J dt'\5B ± (t')\/h~ J dt' cos u>zt'5B x {t'), 

which averages to zero for the low frequency components 
(id < Ldz) Of SBx • 


C. ac Bloch equations 

In the section, we review the master equation in 
the Bloch-Redfield approximation in the rotating frame, 
valid when the Rabi frequency is much faster than the 
RF longitudinal relaxation ratesj^ which is satisfied in 
the parameter regime of interest in this work. These 
equations include the dissipative qubit dynamics due to 
the noise term JB^ and and Sujz in Eq. ( [l6| ). The 
qubit master equation for the pseudospin s = y0jcr|'0)/2, 
where |t/? ) the q ubit state vector, is given by the Bloch 
equation ^ -- ^ , 

ds B 
— = — x s 
dt h 

- D x (s x ~ s x )e x ~ D Y s Y e Y - D z s z e z (24) 

where B is the effective field whose dominant term is 
hojRex (see Eq. ([l6|)), and the relaxation rates are given 
by 


Dx = 

-Llp 

B y = ^7 + 


Dz - + 


Szz{wr) 
2 h 2 

2 h 2 


(25) 


The RF Ti and T 2 relaxation times are given bj @ 


1 _ Syy(u)r) + Szz(^r) 

Ep ~ 2^ 

1 _ D' y + D' z 1 1 

E P ~ 2 


2T ip + TL ’ 


(26) 


where the RF pure dephasing time is given in Eq. (23). 
Unlike the relaxation rates for dc rotation, these rates 
depend on both the Rabi and qubit frequency. 

A plot of the relaxation times (T^,Ti p ,T2 P ) for both 
types of driving as a function of detuning is shown in 
Fig. [6J; and an illustration of the associated relaxation 
processes is shown in Fig. [6i. For optimal parameters, 
Gzz C 1 as shown in Fig. 6k, so that Tip, the decay 
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FIG. 7. (Color online) (a) Illustration of drive amplitude noise 
(solid line ) mediated by virtual transitions (dashed lines) due 
to a combination of ac drive and detuning noise, described by 
second order term Sh ac Eq. Q in the effective Hamiltonian, 
(b) Conceptual illustration of the fluctuations in the drive 
amplitude and associated Rabi frequency 8ujr due to low fre¬ 
quency detuning noise Sc. (c) For (A?,A°) = (30,30) /xeV, 
estimate of the mean square relati ve n oise fluctuations in 
the Rabi frequency, (< t* r /ujr ) 2 [Eq. ([28])] . Blue, solid (gold, 
dashed) curve indicate detuning (coupling) drive. 


time of ac X rotations, is mainly due to noise power in 
the transverse components Sxx oc Gyy, see Eq. (21), 
which scales quadratically with the detu ning drive Rabi 
frequency as mentioned in section MV^ as a draw¬ 
back of detuning drive. This relaxation rate dominate 
the decay near the detuning sweet spot at e* ~ 20 /ieV, 
where T2 P ~ 10 ns. As experimental evidence for this 
relaxation mechanism, we find that for the parameters of 
Ref. Uij X2 P (e*) = 5 ns, consistent with the short relax¬ 
ation times found therein. On the other hand, beyond 
about e = 200 /ieV, decoherence from detuning noise is 
strongly suppressed, with all relaxation times above 1 fi s. 
At this point, phonon dephasing with T p h = 1 /is becomes 
the limiting relaxation mechanism. 

The equilibrium pseudospin s in general depends on 
the response of the noise bath to the driven qubit. Ana¬ 
lytic expressions relating s to the noise power spectrum 
are given in Ref. H] and uni For the numerical simula¬ 
tions in section |VI| we take the RF equilibrium condi¬ 
tion s = — ex,^r corresponding to equal populations 
of qubit states |=b Z) in the lab frame, which is valid in 
the low temperature limit 2k rT <C huiz,R, appropriate 
for typical temperatures T ~ 0.1K of experiments, and 
the regime <C c or < cjz, appropriate near the 

optimal working point. 


Fig. [4]i and Eq. (Cll). An explicit calculation of the 
resulting low frequency decay envelope is given in ap¬ 
pendix |C 3 b| Since we have already tuned the tunnel 
coupling to minimize Scuz = ^eduJz/de (Gzz < 1), the 
quadratic term in Suj 2 z /2ujr = Gzz ^ 2 is negligi¬ 
bly smallThe main threat comes from low frequency 
noise in Sujr, which causes decay of RF X rotations given 
by [cf. Eq. (C19)] 


\sz(t)J 

where 



Wif = exp 


(27) 


°R(t) = \J (&4)lf = > (28) 


the angular brackets (... ) lf denote an average over low 
frequency noise, and the quasistatic detuning noise vari¬ 
ance is defined by [cf. Eq. (C15)] 


r^/t d/ j 

a;(t) = (Se\ = 2 J —S £ (u;) = 


= ^ln 
7r 



(29) 

In the limit uoit <C 1 relevant to the qubit gate times 
t ~ 1 — 10 ns, the logarithmic corrections dominate the 
decay. 

Eq. ( [29]) will be included exactly in the simulations of 
sectionjyij Here, we first give an estimate for the size of 
its effect on the gate infidelity, given approximately by 


1 - Wif (t g ) 



(30) 


where t g oc 1/ljr is the gate time, the effective detun¬ 
ing noise variance, cr* R = a*dujR/dSe is related to the de¬ 
phasing time scale T£(X) = \/2h/cr R , where cr* is given 
in Eq. ©• This infidelity estimate, for the optimal tun¬ 
nel coupling parameters, is plotted in Fig. [7} It shows 
that for large detunings e ~ 200 /ieV, low frequency noise 
causes < 1% infidelity for detuning drive and < 0.5% 
infidelity for coupling drive. We note this estimate is 
most likely an overestimate, since recent experiments in¬ 
dicate that the charge noise spectrum has a milder low 
frequency singularity than 1//. Specifically, Ref. [14] finds 
a spectrum of the form S e (u) oc 1 /cwith /3 = 0.7. 


D. Low frequency dephasing from noise in Rabi 
frequency 


In this section, we present our model for quasistatic 
dephasing due to noise in Sujr and Scuz that is not cap¬ 
tured by the Bloch relaxation rates Eq. (26) of section 
V C| These noise terms cause low frequency fluctuations 
of the instantaneous Rabi frequency Scur + 8uj%/2 ujr 
that is linear and quadratic in 5e, respectively) 8 1 9 1 1Q ] see 


VI. GATE FIDELITIES 

In this section, we compute and optimize the quan¬ 
tum process fidelity of the ac X n (NOT) gate for both 
detuning and coupling driving. We first consider the 
state fidelity of an X n gate, defined as the probability of 
reaching the target state rotated by X n from the initial 
state | —Z). Since the total decay can be factorized into 
the product of the low frequency decay envelop Wu(t) 
Eq. ( |3Q| ) and the exponential decay envelop associated 
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FIG. 8. (Color online) (a) Rabi frequency as a function of detuning e, from Fig. [5j shown here for comparison with infidelity. 

(b) Estimate of the state infidelity of X n rotations (31) based on the gate time, RF Bloch dephasing time T 2p Eq. ( |26| , and 
the dephasing envelop Wu(t n ) Eq. (27). (c and d) Process infidelities for the ac X n gate 1 — F(X n ), where the fidelity F is 
defined above Eq. ( |35| ), as a function of detuning e and drive amplitudes A e and Aa, with colors scaled to Log (1 — F(X 7r )). 
For detuning drive with amplitude A e (c), an optimal point occurs at (e, A e ) — (316,30) peV with F(X 7r ) =s 99.1%. For tunnel 
coupling drive with amplitude Aa (d), an optimal point occurs at (e, Aa) = (400,10) /mV with F(X 7r ) = 99.8%. Both plots in 

(c) and (d) are computed with the static tunnel couplings A? = A° = 30 peV. 


with the RF Bloch relaxation time T 2p ®, the state fidelity 
is given by 


where t n = tt/ujr. Equivalently, we implement Eq. (34) 
on the density matrix as 


F S {X„) = e-^'Wufa) (31) 

where t n = tt/ujr is the X n gate time. The infidelity 
1 — Fs^X^) as a function of detuning for both type of 
driving are plotted in Fig. [8^b), which suggests that fi¬ 
delities exceeding 99% for both types of drive, and that 
coupling drive could yield fidelities near 99.9%. 

We next describe our procedure for numerical simula¬ 
tions of the X n gate in the rotating frame. For comput¬ 
ing the process fidelity, it will be more convenient to use 
density matrix master equation, 


P = p] ~ D 


(32) 


where Hrr is the RF Hamiltonian Eq. (15), without tak¬ 


ing the RWA, D is the dissipator given by 


D = D x {s x - s x )cr x + D Y s Y a Y + D z s z cr z , (33) 

and s(t) as Tr(p(t)<r)/2 is the qubit pseudospin, 
[cf. Eq. (|24|)]. After obtaining the density matrix p(t n ) 
after an X n gate by numerical integration of Eq. ( [32] ) , 
we incorporate phonon-induced singlet triplet dephasing 
described in section ra and low frequency dephasing 
due to Sujr described in section ra which cause the 
pseudospin components to decay as 

Sx{tn) e _t 7 r / Tph Sx(Gr) 

S Y (U t) —^ e~ t 7 V / Tph Wtf(tn)s Y (tn) 

Sz(tn) Wtf(t 7r )s Z (t 7r ) (34) 


p{t ,r) ->• e tv/Tph (s x a x + W li (t n )s Y a Y ) + Wu(t n )s z az ■ 

To find the optimal working point, we compute the 
quantum process fidelity of the X n gate as a func¬ 
tion of detuning and ac drive amplitudes. For a gen¬ 
eral quantum process £, the process fidelity is given by 
F{£) = Tr[x(£)x(£o)]> where x is the process matrix de¬ 
fined b}® 


£(p) = E EmpE n Xmn • (35) 

rn,n 

In our case, the process £ is the X n gate, £{p) is the 
final density matrix after an X n gate computed from 
the simulation procedure described above, for an arbi¬ 
trary initial density matrix p, while the ideal process, 
denoted by £q, is a perfect X n gate in the rotating frame 
£ 0 (p) = cTxpcrx • We follow the procedure for comput¬ 
ing Xmn given in Ref. |46] with a basis set given by 
Em = {1, &x, &z}- 

The X n gate infidelities for detuning and tunnel cou¬ 
pling drive are plotted as a function of detuning (e) and 
drive amplitude (A u ) in Fig. [8^ and d. In large regions of 
the parameter space (e, A w ), detuning drive fidelities ex¬ 
ceed 99% and tunnel coupling drive fidelity reach 99.8%. 
These optimal regions are determined by the competition 
between i) gate speed, which favors small e and large A u , 
ii) decoherence, which favors large e, and iii) nonlinear ef¬ 
fects, which favors small A u . The regions of high fidelity 
occur along the diagonal because as e increases the drive 
amplitude needs to increase to maintain the same gate 
speed, which would otherwise decrease at fixed drive am¬ 
plitude as shown in Fig. [8^i The maximum fidelities are 
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reached at (e, A e ) = (316,30) /leV with F(X n ) = 99.1% 
for detuning drive and (e, A/±) = (400,10) /ieV and 
F(Xk) = 99.8% for tunnel coupling driving. 

As a check on our decoherence model, we performed 
simulations of q ubit dynamics using the three state 
Hamiltonian Eq. (A4), with 1/f detuning noise numeri¬ 
cally generated by the procedure described in appendix 
|C 4| The maximum fidelity for tunnel coupling driving 
computed with these simulations agree with the result 
based on Eq. (l32|) and (1341) . 


VII. DISCUSSIONS AND CONCLUSIONS 


Coppersmith for guidance and stimulating discussions. 
This work was supported by the Intelligence Community 
Postdoctoral Research Fellowship Program. 

Appendix A: Derivation of effective Hamiltonian 

In this appendix, we derive the effective Ha milton ian 
for the qubit, following the standard procedure,^® 0 in¬ 
cluding modifications due to dynamics. We first identify 
the leakage state that is separated from qubit subspace 
by a finite gap at all detunings. To this end, we first diag¬ 
onalize the {| • T), \S •)} subspace, where the tunnel cou¬ 
pling A 2 causes hybridization of the 5(2,1) and T(l,2) 
states. The charge and spin hybridized eigenstates are 


The theoretical framework we have developed for an¬ 
alyzing and optimizing ac gate fidelities in this work 
are quite general, and can be applied to a generic noise 
power spectrum and to other semiconductor quantum dot 
qubits, such as the singlet triplet qubitWhile we con¬ 
sidered silicon, the three-electron double dot hybrid qubit 
can be implemented in other materials such as Germa¬ 
nium or Gallium arsenide.^ Our theory will be applica¬ 
ble in these materials provided that the qubit and noise 
parameters are adjusted appropriately. For example, the 
qubit frequency will be set by the single dot singlet triplet 
splitting in these materials. The phonon relaxation rates 
depends on the character of phonon excitations, attenua¬ 
tion rates, and the lateral dot radius, which is set by the 
transvers e effe ctive mass and gate-defined confinement 
potential ! 41 * 43 * 

While we were focused in this work on how charge 
noise affect qubit fidelities, the decoherence model we 
have developed can be used to turn the question around, 
to probe the charge noise spectrum by measuring qubit 
dynamics P The information thus gained about the na¬ 
ture of the environment causing the charge noise can be 
helpful in developing experimental and fabrication tech¬ 
niques to reduce it. 

In summary, we have systematically analyzed gate fi¬ 
delities of the ac driven quantum dot hybrid qubit, in¬ 
cluding decoherence due to 1// charge noise, determined 
the optimal parameter regime for the tunnel coupling, 
detuning, and ac drive strengths, and showed that gate 
fidelities up to 99.8% can be achieved by driving the tun¬ 
nel coupling. The fidelities computed in this work are 
exponentially sensitive to the 1/f detuning noise param¬ 
eter (c e ), so that we expect even modest reduction in 
charge noise could result in significant further improve¬ 
ments in qubit fidelity. The ac driven single qubit opera¬ 
tions studied in this work are a crucial part of two-qubit 
gate sequences such as the one proposed in Ref. m and 
should thus enable a high fidelity universal gate set for 
quantum dot hybrid qubit. 

acknowlegements -We thank Mark Friesen and S. N. 


(!)j = ( sinr lo cosrjo \ f\■ T)\ 
IjL),/ \cosrjo -sinr/oj \\S ■)) 


where 

f cos77o\ _ 1 | VF 12 + Est ~ e 0 \ / Ai n 

V sin W “ Vv^-^ST + eo) ’ 

770 is the eo dependent mixing angle, and 

E 12 = V / 4Al(eo) + ( J E S T-eo) 2 , 

is the energy splitting, and we define the detuning de¬ 
pendent tunnel coupling A 2 (eo) = A 2 e~^ e °~ EsT ^ €2 . The 
mixing amplitudes (cos 77 , sin 77 ) as functions of detuning 
are plotted in Fig. [9] 

With this change of basis at eo, defined by the following 
transformation 

(10 0 \ 

Urj(e 0 ) = 0 sin 770 cos 770 (A2) 

\0 cos 770 — sin 770 / 

the Hamiltonian becomes 



e (AveV) 


FIG. 9. (Color online) T he m ixing amplitudes cos 77 = (S- 11) 
and sin 77 = (-Tjl) in Eq. (Al) as a function of detuning e. 
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H v (e,A 1 ,A 2 ) = UlH 0 U v = 


( e 0) 

2 

Ai(t) COS 7/0 
y-Ai(t) sin r/o 


Ai(£) cos 7/0 

^ cos 27/0 - A 2 (t) sin 2t] 0 

E ST -e(t ) s * n 2 ^ 0 _ cos 2^ 0 


—Ai(t) sin 7/o \ 

E S T-e(t) sin 2^ 0 _ A 2 (t) COS 27/0 5 

A 2 (t) sin 27/o - cos 27/ 0 / 

(A3) 


where, turning on both detuning and coupling drive, 

e (t) = e o + Cac(t) , Ai(t) = A^ + AJ c (t), A 2 (t) = A 2 (eo) + A^it ). 


Writing = H^ c + H2 C , the dc Hamiltonian is given by 


J*£(e 0 ,A?,A°) 



eo 

2 

A 3 cos 7/0 

-A? sin % 

\ 


A? COS T]o \ 

(^ sr - e° 2 ) 

0 


V 

-A? sin 770 

0 

1 (-Est + E ° 2 ' 

), 


and the ac part is given by 


(A4) 


HUt) 




( -tac/2 Af COS 7/o 

Af c cos 7/0 (e ac /2) cos 2t/ 0 - A^ c sin 2t/ 0 

y—Af c sin t/ 0 -(e ac /2) sin 2t/ 0 - A§ c cos2t/ 0 


—Af c sin 7/o 

-(e ac /2) sin 27/ 0 - A^ c cos 2 t/ 0 
Ag c sin 27/o - (e ac /2) cos 2 t/ 0 


Detuning noise can be included simply by taking 
e ac e ac + Se. So far, no approximations have been 
taken. Note that has different dependences o n e a c 
and eo, due to the eo-dependent transformation Eq. (A2). 

Next, we perform a canonical transformation to de¬ 
rive the effective Hamiltonian for the ac driven qubit. 
This transformation is a perturbative expansion orga¬ 
nized as follows. We first separate H v into a diago¬ 
nal H d = diag (H) and an off-diagonal = H — H d 
part, which is considered to be a perturbation, so that 
H v = H d + The off-diagonal part is further sep¬ 

arated into tt= Vq + 14, where Vq is off-diagonal 
within qubit subspace, Vl is off-diagonal between qubit 
and leakage space. Vq is justified as a perturbation even 
for large Ai because in the far detuned regime cos t/o«1 
(see Fig. [9]). We then perform the canonical transforma¬ 
tion 


I ip) = U s \i>), U s = e s , 

where S is anti-Hermitian (S^ = —S) and a is purely 
leakage term. The transformed Hamiltonian, up to 
0(S 3 ) terms, is given by 

H = UgHU s - ihU ] s d t Us 

• ifi • 1 if) - 

= H — ihS + [S,H + ^S] + -[5, [S, H + jS]} 

1 ih 

= H d + V Q + [Vl, 5] + ~[[H d , S], S] + -[S, S} 

+ V L + [H d + V Q ,S] + ihS + |[[V S],S] (A5) 


The lea kage terms, given in the last line of 
Eq. (A5), are eliminated perturbatively in the pa¬ 


rameters a n = /AE n i and /3 n = fvuoz/AE n i, where 

A E n i = E][ — E^ E ^ t are the diagonal elements of H v , 
n(l) labels states in the qubit (leakage) subspace. For 
static Hamiltonians, a n is sufficient to parametrize the 
perturbation theory. In the dynamic case, the additional 
expansion parameter (3 comes from the dynamics de¬ 
scribed by S. We will denote the kth order of the pertur¬ 
bative series as the term of power with k = p + q. 

The decoupling procedure is iterative® S = Y^kLi 
is a perturbative series with each term of order &, chosen 
to eliminate leakage terms of order a^ _1 , and produces 
corrections in the effective Hamiltonian as well as leakage 
terms of order k. 

In the leading order, S = is chosen to satisfy 

[H di SW] = -V L , S^ = -VA, 


The leakage term to the next order, generated by is 
given by 


v« = [v Q ,s {1) }+ihd t s (1) , 


which can be eliminated with ,S' (2 * = —V^/A E n i, re- 
suiting in terms of order <a m <a n and <a m/ d n in the effective 
Hamiltonian. 

We will keep the perturbation theory to leading order, 
where the nonzero matrix elements of = S dc -\-S ac (t ), 
S dc being the static part and S ac the time dependent 
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part, are given by 


The effective Hamiltonian is then given by 


< 2 | S dc | 0 ) = 


2 A® sin r/o 


<2| S ac (t) |0) = 


A“ c sin r/o 


(2\S ac (t)\l) = - 


A£'o2 

(e/ 2 ) sin 2r/o — A !) 0 cos 2r]o 


A E, 


02 


A E 


12 


(A 6 ) 


h = P ri HP v = H d + V Q + l -[V L , S&] , (A7) 

where P p = \rji) (rji \ + I 772 ) (7721 is the projection operator 
onto the qubit subpsace, defined as the lowest two energy 
eigenstates 1 770 , 1 ) of H r] at t = 0. The matrix elements 
are given by 


k ( t ) = _ 2 A^(t)sin 2 7/0 _ 

<Unl ' 2 cos 2t]q(Est ~ e(£)) + 2 A 2 (t)sin 2 r / 0 + E S t + e(t) 

, /.N , 1 Am/^ , cos Vo(Est — e (t)) + 2A 2 (t)sin77o 

Aoi W = ftio = 7 A iW cos 770 +-- 

4 V cos 2770 (# 5T 


^11 W = 7 ( 2 ^ 5 t - 0 

4 V cos 2t]o(E S t 


4cos tjoEst + 4A 2 (t)sin?7o 


e(t)) + 2A 2 (£)sin2?7o cos2t]o(E S t 
A 2 (t) + (i^sT - e(t)) 2 


5(t)) + 2A 2 (£)sin2?7o 


e(t)) + 2A 2 (t)sin 2ryo + 

+ cos 2rjo(e(t) - E S t ) - 2A 2 (t)sin 2r/ 0 ) , (A8) 


r 


We note here that the difference in the eo and e ac depen¬ 
dence comes from the implicit eo dependence in 770 , and 
in particular the difference in the respective derivatives 
is given by 


dh dh 
de 0 de ac 


dr] d 
de 0 dr\ ’ 


(A9) 


where 77 is the mixing angle defined in Eq. ( Al| ). This 
difference is important at low detunings, where r] has 
strong detuning dependence, as shown in Fig. [9] 

The static effective Hamiltonian that follows by taking 
e(t) eo is given by 


H ( 


dc 


e 0 /2 - j A i s 7 2 7 


Ai cos 770 
and the qubit basis is 

I Qn) — (1 - SdJU^qn) 


Aicosr?o \ 

E$t ~ E12/2J 


(A 10 ) 


(All) 


where q n ) labels the basis states {| • S), |5 •), | • T)}. 


Appendix B: longitudinal driving 


The longitudinal driving field in Eq. (10) is given by 
( dBz 


B 


5e 


d 2 B z 


acz(t) = A u COSUzt 


(El) 


Consider first the dominant contribution in Eq. (Bl) 


proportional to dBz/du , plotted in Fig. 10 r as a func¬ 
tion of e for both types of driving, This term can 
cause modulations and changes in Rabi frequency when 
ujr ~ uoz- However, near the optimal working point of 
this qubit, ujr c oz, so that these effects are strongly 
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FIG. 10. (Color online) As a function of detuning e: (a) The 
longitudinal drive coefficient for detuning and tunnel coupling 
drive, dBz/dg. (b) The coefficient of the noise fluctuation 
the longitudinal drive amplitude d 2 Bz /dedu for both types 
of driving. 


suppressed. For detuning driving, this term, given by 
dBz/de = dEio/de, was already minimized in 


III 


and 

is essentially zero, As shown in Fig. [lOja). For cou pling 
driving, dBz/dA ac ^0.1. The second term in Eq. (Bl) 


represent noise in the drive amplitude that has the same 
origin as Sujr, the noise in the Rabi frequency. In prin¬ 
ciple, this noise term should be added to SBz , which 
could cause Ti P relaxation. However, it is strongly sup¬ 
pressed because: i) Due to the prefactor cosu/^t, the 
noise spectrum is shifted similarly to SBx • In particu¬ 
lar, the noise power that would contribute to Ti p occurs 
at S € (ujz ± wr), which is reduced by an order of mag¬ 
nitude from S € (ujr ), and ii) The drive noise coefficient 
d 2 Bz/dud(5e), plotted in Fig. 10 3 , is of order 10 -3 . 





























14 


Appendix C: 1/f noise 

In this section, we summarize properties of 1// noise. 


1. Noise correlations 

A parameter V with Gaussian noise SV is completely 
characterized by its autocorrelation function, defined by 

/ OO 7 

7^ eW-^Sviw) , 

-oo 

where Sy(oo) = (|(W(c<;) 2 |) is the noise power density. 
The variance isgiven by 

4 = (5V\t)) = S v (t = t') = l°° E S v (w) , (Cl) 

and the correlation time is defined by 

i r°° 

t v = — dt'S(t f ). (C2) 

°V Jo 

For a noise source with exponential correlations, 
S(t) =cTye _t / rv . For 1/f noise, a high (c Jh) and low 
(uji) frequency cutoff is necessary to make ry and cry fi¬ 
nite. The low frequency cutoff sets the correlation time, 
Ty = 27 l/uJi. 


2. Experimental fit of detuning noise strength 

We fit the single parameter c e in the noise spectrum 
using the T\ relaxation time for decay from |1) -4 |0), 
measured to be 7 ns at the sweet spot in Ref. [211 A sim¬ 
ilar time scale was measured in Ref. Ena which attributes 
it to charge noise. We thus fit the detuning noise strength 
to the T\ relaxation time using the Bloch formula 

1 _ Sxx{^z) _ GxxS e (io z ) 

Ti “ 2 h 2 ~ 2 h 2 

from which we find 


I 2 h 2 id Z (e*) 

Gxx(e*)7i(e*) 


= 2.38 /ie v , 


(C3) 


where the detuning sweet spot for parameters in 2 D is lo¬ 
cated at e* ~ 20 /ieV. 


3. Pure dephasing rates due to 1/f noise 

In this section, we review the relevant formulae for 
pure dephasing rates for 1/f noise. Consider two 
quantum states denoted generically as | ± ), with an 
energy splitting E = E+ — E_ and noise fluctuations 
SE = 5E+ — 5E -, which for definiteness we assume in 


the section to come from detuning noise Se. Pure de¬ 
phasing refers to the decay in the off-diagonal element 

of the density matrix p H = ( + \p\ — ) = e~ %< ^ /2, 

where </> = (E+ — E_)t/h is the relative phase and 


5(/>(t)= [ dt'5E(t')/h , (C4) 

Jo 

is the accumulated phase due to noise fluctuations, which 
typically causes decay when averaged over noise realiza¬ 
tions. We keep to quadratic order in detuning noise, so 
the energy fluctuation is given by 


5E(t) = aiSe(t) + y 5e(t) 2 , 
( dE d 2 E 

~ \d{6e)'~d{8¥ 

The decay envelope is given by 



W(t) = ( e -^ (t) ) = 


exp 




(C5) 


(C6) 


where T^\t) and T^\t) are due to noise averaging over 
the linear (Se) and quadratic (5e 2 ) term, respectively. 
The linear term can be expressed in a well known, closed 
form 


<1 \ t ) = lf o J dhdhiSE^SEih)) 

= - 7 r f [ dtidt2(5e(ti)6e(t 2 )) 
z Jo Jo 

aH 2 duo 2 

= 2 wL 00 to S * u)waD M/2) ’ 


(C7) 


(C8) 


where sinc(x) = sin(x)/x. Note that due to the sine func¬ 
tion, the integral is dominated by the quasistatic part of 
the spectrum (u < 1/t). Physically, this stems from the 
fact that the phase accumulation 5(f) coming from noise 
at frequenci es h igher than 1/t tends to time-average to 
zero in Eq. (C4). 

The decay exponent due to quadratic fluctuations can¬ 
not in general be expressed in closed form, but it can be 
expressed in terms of a functional determinant 


r^ 2) = - In [det(l + ia 2 s/h)\ 1 / 2 
= i tr ln(l + ia 2 S/h) 

= ^itr[(-ia 2 S/S)1 (C9) 

n 


where the multiplicaton of noise correlation functions de¬ 
notes integration over functional kernals 


tr S = 


j dt'S(t/t '), 
tr S 2 = J dt'dt"S(t\t")S(t'\t'). 
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The series expression in Eq. (C9) was previously derived 
in Ref. [T9] , and it can be simply generalized to include 
dynamical decoupling sequences such as spin echo or 
CPMCP. 

At a fixed gate time t g , the summation Eq. ( |C9| ) can 
be evaluated for high and low frequency ranges relative 


to 1/tj 19 . In contrast to the linear term Eq. (C8), both 


high and low frequency ranges contribute at all times, 
but the decay is dominated at short (long) times by the 
low (high) frequency contribution. Here, the short (long) 
time regime is defined relative to the time scale 


T2 = 


h 

a 2 c 2 e 


(CIO) 


b. Dephasing due to low frequency noise 

Next, we consider t he q nasistatic contribut ions to de¬ 
phasing for linear Eq. (|C8|) and quadratic Eq. (C9) terms 


in t urn. Due to 1// singularity, the decay exponent in 
Eq. (C8) is dominated by low frequency noise at uj 1/t, 
so that we can take sinc(cj£/2) ~ 1 in the integrand. This 
yields a Gaussian-like decay given by 




(C14) 


(C15) 


For the optimal parameters considered in this work, 
we will be in the short time regime t <C r 2 , and the 
quadratic noise terms are strongly suppressed. However, 
these higher order effects may be important in the exper¬ 
iments that operate away from these optimal parameters, 
for example, in Ref. [21] and 24. 


a. Linear and quadratic noise couplings 


When the low frequency term (C14) dominates the decay 
(C6), the dephasing time scale is defined by 1 = T^(T 2 *), 


and is given approximately by 19 ^ T 2 « r\ / y/ln (r m /t\ ), 
where t\ = \f2h/a\c e . If we define a variance by the time 
scale T 2 * = %/2h/aia*, then 


= c e \J (1/7r) ln(r m /ri) ~ c e \/(l/7r) In (c e T m /h) 

= 5.5 /jcV, (C16) 


In this section, we determine the linear (ai) and 
quadratic (a 2 ) couplings to detuning noise [Eq. ( |C5| )] 
coming from the leading order expansion in Se of the 
noise field given in Eq. <§• These coefficients can be 
used in the formulae given in appendix |C 3b| and |C 3 c 


to compute dephasing envelops for dc Z rotations in the 
lab frame and ac X rotations in the rotating frame as 
described 

In the lab frame, the fluctuations in the qubit fre¬ 
quency (dc Z rotation frequency) is given to quadratic 
order by 


/- SB Z 

5E = y (B z + SB z f + SB\-B Z =5B Z + ^L, 


For the ac gates times of t g ~ 1 ns, the cr' e (t g ) = 5.9 /xeV, 
which is approximately equal to cr*. 

The dephasing exponent due to the quadratic term for 
a 1// noise spectrum Eq. ( |C9| ) is evaluated in detail in 

Ref. [T9[ where it is shown that the low frequency contri- 
( 2 ) 

bution to Ty is given by 

d 2) W = M 1 + ia 2 taf{t)/h) . (C17) 

This is the dominate contribution for times up to 
t ~ r 2 j ln(r m /r 2 ). For short tim es t <C r 2 , it gives a 
Gaussian-like decay similar to Eq. ( |C15 ). 

To summarize, the low frequency (qu asist atic) contri¬ 
butions to the total decay envelope Eq. ( |C6| ) is given by 


where Bz = hwz , and we define Sf s 5edf/d{5e) for a 
generic function of detuning f(Se). It follows that 


ai(Z) 


dB z 
d(Se) ’ 


( dB x /d(6e )) 2 
a 2[^) = -^- • 


(Cll) 


For ac X-rotations, in the RF and RWA, the fluctuations 
of the Rabi frequency is given by 


SE 

h 


f (wr + 6 wr) 2 + <5w| 


■ WR 


5u r + 


2 ur 


(C12) 


so that 


ai(X) = h 


Sljr 

Se 


a 2 (X) 


h 5ui 2 z 
u>r Se 2 


A u d 2 Bx 

2 d{6e)du’ 

= ( dBz/djSe )) 2 

hu) R 


(C13) 


W u (t) = 


exp [—(aitcr'(t)/\/2h) 2 ] 


■\/l + ia^o^jpjh 
For f -C T 2 - the decay envelop is given by 
|Wif(t)| = exp [-( 7 c/,(t)t) 2 ] , 
where the decay rate is given by 




V2 


+ 


a 2 (j' 2 {t) 


(C18) 


(C19) 


(C20) 


When the high frequency cutoff is below the releva nt 
gate speeds, uoy <C 27r/t, the decay envelop Eq. (C18) is 
equivalent to the one computed from a Gaussian aver¬ 
age over static noise with the variance given by the total 
integrated noise power, 


where u = e ac or A ac 


<r e = ( Se 2 (t )) = c e \j (l/7r) ln(u h /uji ), (C21) 
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X rotations Z rotations 



0 50 100 150 200 0 50 100 150 200 

e (jueV) e (£/eV) 


FIG. 11. (Color online) The exponential decay time scale T 2 
as a function of detuning: (a) For ac X rotations, with the 
tunnel couplings(A5, A°, r) m (30,30) //eV, for detuning (e) 
and tunnel coupling (A) drive, (b) For dc Z rotations, with 
the tunnel couplings (A?, A°,r) = (30,30) /xeV used in this 
work and (A?, A°) = (10,30) /xeV used in Ref. [24. 


which can be significantly larger than a*. For example, 
for uj h = 100 GHz, a e = 9 /ieV. However, as discussed 
in appendix | C 2[ the noise spectrum is finite at GHz fre¬ 
quencies, so the limit ujh 27t/£ is not satisfied. 


c. Dephasing due to high frequency noise 



t (ns) 



FIG. 12. (Color online) (a) 1/f detuning noise generated nu¬ 
merically. (b) The probability of occupying |1) as a function 
of time (dots), computed by numerical integration of the mas¬ 
ter equation with detuning noise generated as shown in (a). 
(Solid line) Fit to an exponential decay with T\ — l ns. 


t 2 (Z) oc Gxx Wr, [cf. Eq. (Cll)] scales quadratically 
with the detuning driven Rabi frequency, which was in¬ 
creased in going from the tunnel coupling Ai = 10 /ieV 
(Ref. U1J to Ai = 30 /ieV (optimal) in this work. 


The high fre quen cy component of the quadratic noise 
term 5e 2 (Eq. ( |C9| )) cause s an exponential decay e _t / T2 , 
where r 2 is giv en i n Eq. (CIO), so that the total decay 
envelop in Eq. (|C6|) is given by 


(C22) 


W(t) = e-^Wuit). 


We plot the exponential decay time scale 72 (X) for ac 
X rotations in Fig. [Tl] This decay time is > 10 /is near 
the optimal working point (e > 200/ieV), which, for the 
typical ac gate times considered in this work of < 10 ns, 
causes infidelities of < 0.1%, and is thus negligible. 

The exponential decay time scale 72 (Z) for dc Z ro¬ 
tations are plotted in Fig. mi for both optimal tunnel 
couplings considered in this work and that of Ref. l24l 
Although this decay time 72 (Z) = 1 /is near e ~ 200 fieV 
is significantly shorter than 72 (X) rotations, due to the 
short Z gate periods ^ 1 ns, Z the gate infidelity, 
which can be estimated as t g (Z)/r 2 (Z) ~ 0.01%, is still 
very small, as noted in section m Note that this de¬ 
cay time is actually shorter for the optimal tunnel cou¬ 
plings in this work then that of Refj24] This is because 


4. Simulations with numerically generated 1/f 
noise 

In this appendix, we describe our simulations of qu bit 
dynamics using the three state Hamiltonian Eq. ( |A4| ), in 
the presence of numerically generated 1/f detuning noise, 
following the procedure described in 50 . We first generate 
a white noise time series /(£), whose Fourier component 

f u = J dt'e^'fit') 

is Gaussian distributed with a uj independent variance, 
which we set equal to the detuning noise variance 

S f ^) = (\U 2 )=aj. 

We then multiply the stochastic Fourier components by 
1 / y/w, and construct a 1// noise time series 

< c23) 
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that has a noise spectrum 





UJ 


which gives the desired noise spectrum if we set white 
noise variance equal to the noise strength, c e = <jf. To 
check that this is the correct noise strength, we simulate 
T\ relaxation by initiallizing the qubit in the state |1) and 
computing the probability to remain in |1) as a function 
time. Fitting this probability to an exponential decay 


Pi = e -t / Tl + (1 — e~ t / Tl )Pif yields the relaxation time 
T\ — 7 ns, consistent with the experimental result used 
to determine c e , see Fig. [I2| 

In our simulations, we find it convenient to use the 
discrete cosine transform, where all Fourier components 
are real from the outset (instead of the discrete Fourier 
transform). In computing the qubit fidelity, the ideal 
gate is defined as the X n ac gate operation without any 
noise, see Eq. (35). We find that the qubit fidelity con¬ 


verge after averaging the solution to the density matrix 
equation of motion over 20 realization of the Se(t) time 
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